Detection of myocardial ischemia by magnetocardiogram 
using 64-channel SQUID system 

M. Sato 1 , Y. Terada 2 , T. Mitsui 2 , T. Miyashita 3 , A.Kandori 3 , and K. Tsukada 3 

department of Thoracic Surgery, Fujita Health University School of Medicine, Aichi, Japan; 2 Department of 
Surgery, Institute of Clinical Medicine, University ofTsukuba, Ibaraki, Japan; 3 Central Research Laboratory, 

Hitachi Ltd., Tokyo, Japan 


1 Introduction 

Ischemic heart disease (IHD) is one of the most 
critical diseases in view of its serious mortality and 
frequency of occurrence. Especially in Japan, aging 
population and recent changes in life styles, such as 
increasing intake of animal fat, are expected to bring 
about a more frequent occurrence of atherosclerotic 
lesions in coronary arteries, resulting in an increased 
IHD morbidity. Furthermore, we can now see more 
IHD in younger population. Therefore, accurate and 
quick diagnosis is important. For early detection of 
IHD, a non-invasive and convenient diagnostic tool 
with reasonably high accuracy is needed. From this 
point of view, resting standard electrocardiography 
(ECG) and echocardiography, which are standard 
non-invasive diagnostic tools, are not satisfactory in 
the diagnosis of IHD. 

The electrical current generated by the heart creates 
not only electric potential but also a magnetic field. 
We have observed electrophysiological phenomena 
of the heart by measuring tangential components of 
magnetocardiograms (MCGs) using a multi-channel 
superconducting quantum interference device 
(SQUID) system [1]. Using the tangential 
components, we have obtained a two-dimensional 
current projection pattern on the observing plane, 
because these components show a peak pattern just 
above the current source. We have developed vector 
arrow maps of the current from tangential MCGs 
and have observed their time courses in normal 
subjects and various patients with cardiac diseases 
[2,3,4]. We have also developed a visualization 
method, iso-integral mapping, that is used for two- 
dimensional visualization of total current 
distribution during the depolarization and 
repolarization of the ventricle [5]. 

With our multi-channel SQUID system and these 
two visualization methods, we thought we may be 
able to use MCGs as an absolutely non-invasive and 
convenient diagnostic tool for IHD. In this study, 
we tried to visualize myocardial ischemia by 
measuring MCGs using a 64-channel SQUID 
system, and compared the effectiveness of this 


approach to resting standard 12-lead ECG and 
echocardiography. 

2 Methods 

2.1 Subjects 

MCGs were obtained from 25 normal subjects and 
25 patients with IHD. Informed consent was 
obtained from all subjects and patients before 
participation in the study. 

1. Twenty five normal male and female subjects 
aged 21 to 44 years (mean±SD: 30.4±5.8) were 
studied. Their mean height was 1.71 ±0.05 m, 
mean weight was 68.16 ±9.93 kg, and mean 
body mass index (BMI) was 23.17±3.10 kg/m 2 . 
They had no history of cardiovascular disease, 
abnormal physical findings, or signs of 
abnormalities in their standard ECG. 

2. Twenty five IHD patients aged 44 to 74 years 
(mean±SD: 63.2±7.7) were studied. Their mean 
height was 1.60±0.05 m, mean weight was 62.2± 
8.7 kg, and mean body mass index (BMI) was 
24.15±3.37 kg/m 2 . All of them were treated in 
Tsukuba university hospital from December 
1998 to September 1999. They were examined 
using resting standard 12-lead ECG, 
echocardiography, and coronary arteriography 
(CAG). Significant (75% or more) stenosis of 
the left anterior descending artery (LAD) and/or 
left main trunc (LMT) was found in every 
patient by means of CAG. Significant stenosis 
of the left circumflex artery and/or right 
coronary artery existed in 20 patients. The 
resting ECG examination showed ischemic 
changes, such as an abnormal Q wave and ST-T 
change, in 18 patients, but there was no 
ischemic change in seven patients. 
Echocardiography revealed abnormal left 
ventricular wall motion, asynergy, in 16 patients. 



2.2 MCG measurement 

We recorded MCGs using the 64-channel SQUID 
system developed by our group [6] and installed in 
Tsukuba University Hospital. The MCG 
measurement was done inside a magnetically 
shielded room (MSR). The system had 64-channel 
magnetic sensors inside the dewar, and each 
magnetic sensor probe incorporated a first-order 
gradiometer. A probe was placed at each of the 64 
sites in an 8x8 matrix with 25-mm intervals, giving 
a measurement area of 175x175 mm. This area 
could cover the entire heart of an adult. 

Subjects with dresses on lay on a bed inside the 
MSR, all metallic objects were removed. The dewar 
containing SQUID sensors was placed as close as 
possible to the anterior chest without touching it. We 
set each sensor to its exact position by using a laser 
pointer. We measured the normal (Bz) component of 
the magnetic field vertical to the torso for 30 
seconds. Simultaneously, we recorded the ECG lead 
II. The signals could be viewed on a real-time 
display, and the data was stored and analyzed by a 
personal computer. 

2.3 MCG analysis (Vector arrow map) 

The obtained MCG signals were averaged by the 
personal computer. To investigate the location and 
the direction of the current source at each point of an 
MCG waveform, we constructed a combined 
isomagnetic and vector arrow map from spatial 
derivatives of the normal (Bz) component in the 
tangential direction[7]. The imaging current was 
defined by rotating the magnetic field arrow 90 
degrees counterclockwise according to Biot-Savart’s 
law. The peaks of the magnetic fields appeared to 
coincide with the most active regions in the heart, as 
deduced from the tangential MCG characteristics. 
From these peaks, the location of the current source 
was easily recognized. 

2.4 MCG analysis (Iso-integral map) 

Iso-integral maps for two time intervals, the 
depolarization process (QRS) and repolarization 
process (ST-T) which were defined from the MCG 
waveforms, were constructed. The QRS iso-integral 
maps and ST-T iso-integral maps were obtained by 
calculating, for each point, the algebraic sum of all 
magnitudes from QRS onset to the J point and from 
the J point to T wave offset, respectively. Maps were 
drawn by integrating the magnitude at each 
measurement point. 


We obtained MCG data from all the normal subjects 
and IHD patients, reconstructed vector arrow maps 
to show the current distributions for each time, and 
also reconstructed QRS and ST-T iso-integral maps 
to visualize the total current image for each time 
interval. 

3 Results 

3.1 Vector arrow and iso-integral maps from 
normal subjects 

The distributions of current arrows for all subjects 
during QRS and ST-T were similar. Upon onset of 
the QRS wave, current arrows were directed to the 
lower right at the intraventricular septum. These 
arrows changed direction anticlockwise and were 
directed to the lower left at the left ventricle. Then 
arrows were directed gradually to the upper right at 
the basal portion of the heart. During ST-T, current 
arrows were directed to the lower left consistently. 
These results were the same as those of our previous 
study [1]. 

An example of the iso-integral maps from normal 
subjects is shown in Fig. 1. As for all healthy 
subjects, these iso-integral maps show a similar 
pattern for both the depolarization and the 
repolarization process. Both maps show a similar 
oval shape with the maximum value located at the 
position of the left ventricle. The maximum integral 
value of repolarization (j ST-Tmax) was higher than 
that of depolarization (j QRSmax) in each case; 
mean I QRSmax was 15963.86±7763.24 pTms/cm, 
mean J ST-Tmax was 22456.94±9620.30 pTms/cm, 
and the mean j ST-Tmax / J QRSmax ratio was 
1.49±0.37. 
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Figure 1: Iso-integral maps of a normal subject, a) 
Iso-integral map during QRS. b) Iso-integral map 
during ST-T. 



























3.2 Vector arrow and iso-integral maps of 
IHD patients 


whose echocardiographic results did not show 
ventricular asynergy. 


In 22 out of 25 IHD patients (88%), we found 
abnormal mapping patterns during the repolarization 
phase. The abnormalities were classified into three 
types: a rightward shift of current arrows in the 
vector arrow map (n=13, Fig. 2), a multidipole 
pattern in the vector arrow and iso-integral maps 
(n=10, Fig. 3), and a decrease of the maximum 
integral value in the iso-integral map (n=15, Fig.4). 



Figure 2: a) A case showing the rightward shift of 
current arrows, a) Vector arrow map with a 14 msec 
interval, b) Iso-integral map during ST-T. 



Figure 3: A case showing the multidipole pattern 
during ST-T. a) Vector arrow map with a 14 msec 
interval, b) Iso-integral map. 


Mean j QRSmax was 10957.7±5143.7 pTms/cm, 
mean J ST-Tmax was 11160.4±6364.1 pTms/cm, 
and the mean j ST-Tmax / \ QRSmax ratio was 
1.07±0.39. The mean j ST-Tmax and the j ST-Tmax 
/ j QRSmax ratio were significantly lower than those 
of normal subjects (pO.OOOl, p=0.0014, 
respectively by the Mann-Whitney U-test). 

We found some of the MCG abnormalities described 
above in four (57.1%) of seven IHD patients despite 
normal findings for these patients in the resting 12- 
lead ECG examination. We also found MCG 
abnormalities in seven (77.8%) of nine patients 



Figure 4: A case showing the decreased maximum 
integral value during ST-T. a) Iso-integral map 
during QRS. b) Iso-integral map during ST-T. The 
j ST-Tmax / 1 QRSmax ratio was 0.909. 

4 Discussion 

In this study, we found MCG abnormalities during 
the repolarization phase in 88% of patients with 
severe coronary lesions. The abnormalities of the 
vector arrow maps and iso-integral maps could be 
classified into three patterns: a rightward shift of 
current arrows, a multidipole pattern, and a 
decreased maximum time-integral value. 

The first two types, the rightward shift and 
multidipole pattern, were previously reported to be 
the expression of various heart diseases such as left 
ventricular overloading and old myocardial 
infarction[8,9]. Although why these abnormalities 
appear in IHD patients is unclear, they are thought 
to be caused by the inhomogeneous repolarization 
currents. 

In a normal myocardium, the current flow is almost 
uniform during repolarization. When there was an 
ischemic part, the current flow seemed to become 
random. Subsequently, the strength of the magnetic 
field would become weaker. Such a phenomenon 
would decrease the maximum integral value in the 
iso-integral map. 

In this study, almost all patients had multiple severe 
ischemic lesions (including myocardial infarction) 
with three dimensional extension, indicating a need 
for surgical treatment. Because the visualization 
methods we used were based on projection to a two 
dimensional plane, we did not investigate whether 
the MCG could clarify the localization and degree of 
ischemia. 

On the other hand, this study suggested that MCG 
had higher sensitivity for detecting myocardial 
ischemia than resting 12-lead ECG or 
echocardiography. Although further study is needed 

























































































































including MCG measurements of patients with 
single coronary lesions, a magnetocardiographic 
comparison between lesions with and without 
infarction, and determination of the detection rate 
for silent myocardial ischemia[10] and syndrome 
X[11], MCG may have a potential as a usefull tool 
for detecting myocardial ischemia in clinical testing 
especially in screening for IHD. 
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